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^ ■ The study of a self-consistent system of nonlinear spinor and Bianchi type I gravitational 

O . fields in presence of a viscous fluid and A term with the spinor field nonlinearity being some 



(N 



(N 
(N 



o 



X 



arbitrary functions of the invariants I an J constructed from bilinear spinor forms S and P, 



generates a multi -parametric system of ordinary differential equations Lli,i2J. A qualitative 
analysis of the system in question has been thoroughly carried out. A complete qualitative 
classification of the mode of evolution of the universe given by the corresponding dynamic 
system has been illustrated. 
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^ ; I. INTRODUCTION 
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Though the investigation of relativistic cosmological models usually has the energy momentum 
tensor of matter generated by a perfect fluid, to consider more realistic models one must take 
into account the viscosity mechanisms [3^ iJ, H, @, 0]- On the other hand in the recent years 
O ■ anisotropic cosmological models with nonlinear spinor field have been extensively studied due to 
the facts that (i) introduction of nonlinear spinor field into the system leads to the isotropization 
of initially anisotropic universe fsLH.floll: (ii) spinor field nonlinearity in some cases can give rise 
■ to singularity-free solutions \l(X \Vu\ : (iii) spinor field can be considered as one of the possible 
9^. , candidate to explain the late time acceleration of the universe [ESQ. 



Given the importance of viscous fluid and spinor field to model a realistic universe, recently 
we have considered a self-consistent system of nonlinear spinor field and a gravitational field 
described by a Bianchi type I (BI) cosmological model filled with viscous fluid Ini,!] we 

have thoroughly studied the corresponding field equations. Exact solutions of the field equations 
were given in terms of volume scale of the BI space-time T. Multi-parametric system of equations 
for volume scale T, energy density of the viscous fluid e was solved for some special choice of 
bulk and share viscosity. Given the richness of the system mentioned above in this report we study 
it qualitatively for more general cases. In the absence of viscosity the system allows integrals of 
motion, whereas it is either impossible obtain, or there is no first integral at all in general when 
the viscosity is taken into account. As a result the study of possible modes of evolutions becomes 
very difficult. Undoubtedly, the result of the investigation should be presented in the form of 
numerical values and at the same time can not be reduced to a representation of some causally 
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found examples of numerical solutions. Clearly, it should be a classification of modes of the 
evolution in the parametric space. Actually, it is the task of qualitative analysis. 



II. BASIC EQUATIONS 

We consider a self-consistent system of nonlinear spinor and Bianchi type-I (BI) gravitational 
fields filled with a viscous fluid in presence of a cosmological term. As it was shown in the 
components of the spinor field and metric functions can be expressed in terms of volume scale T 
of the BI space-time. So one needs to find the function T, explicitly. Corresponding equation can 
be derived from Einstein equations and Bianchi identity [a detailed description of this procedure 
can be found in [li 0]]. For convenience, we also define the generalized Hubble constant. The 
system then reads: 



T 

H 



3Ht, 
1 
2 



K fm ^ — 2) 



{3^H -co)- {3H^ -Ke-A) + 
e = 3H{3^H-a})+4r]{3H^-K£-A)-4riK 
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T 



(2.1a) 
(2.1b) 

(2.1c) 



Here k is the Einstein's gravitational constant, A is the cosmological constant, A is the self- 
coupling constant, m is the spinor mass and n is the power of nonlinearity of the spinor field 
(here we consider only power law nonlinearity). In (12.11) 7] and £, are the bulk and shear viscosity, 
respectively and they are both positively definite, i.e., 

7] > 0, <^ > 0. (2.2) 

They may be either constant or function of time or energy. We consider the case when 

T7=Ae«, <^=5e^, (2.3) 
with A and B being some positive quantities. For p we set as in perfect fluid, 

p=Ce, Ce(0,l]. (2.4) 

Note that in this case C 7^ 0, since for dust pressure, hence temperature is zero, that results in 
vanishing viscosity. Note that a system in absence of spinor field has been studied in [Il4l.ll5n. In 
that case the corresponding system is analogical to the one given in (|2.1I) without the third terms 
in (llJbl) and (IZTcl) . 



in. QUALITATIVE ANALYSIS 

Research on the behavior of the dynamic system given by a system of ordinary differential 
equations implies the survey of all possible scenarios of development for different values of the 
problem parameters. It is necessary to understand at least how the process of evolution comes to 
an end if it does so at infinitively large time for a given set of initial conditions which can be given 
anywhere. 

So, under the specific behavior of the system we understand the phase portrait of the system, 
i.e., the family of integral curves, covering the total phase space. It is easy to imagine as far as 
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any point of the space can be declared as the initial one and at least one integral curve will pass 
through it (or it will be fixed point). 

Certainly, it is difficult to imagine such a set of curves. In many cases, close (and not only) 
curves transform into each other at some diffeomorphism of space. These curves are known as 
topologically equivalent. The differences between them are not very important for our study. 
They all behave in the same manner. This relation - "the relation of equivalence" - divides the 
family of curves into the classes of equivalence. For graphical demonstration it will be convenient 
to present at least one representative of each class. 

The change of the value of problem parameters not always results in significant change of the 
phase portrait. Repeating this method, we say that one family of integral curves (covering the total 
space) for the given set of parameters is equivalent to the other for another set of parameters, if 
there exists a diffeomorphism of space transforming the first family into the second. It is clear 
that there occurs the division into the classes of equivalence, and we are not very interested in 
differences between equivalent families. We argue that the corresponding changes in parameters 
do not alter anything on principle. So it is sufficient to demonstrate only one phase portrait for a 
given set of parameters underlining the features of the given class. 

However, for some critical relations between the parameters there occurs significant changes. 
These are the boundary relations of parameters, dividing, as usual, parameter space into regions 
of similar behavior. Thus accomplishes the qualitative classification of the mode of evolution of 
dynamic system. Now, giving the concrete value of parameters, we can define which region of 
parameters they correspond to, thus define the type of behavior. Moreover, given the specific 
initial conditions, we can answer the question to which region of phase space the evolution of the 
system lead in time. 

In our cosmological model, numerical parameters A, a, B, /3 are related to the viscosity, while 
X and A are the (self)-coupling and cosmological constants. 

Initially, we consider the system of Einstein and Dirac equations. Solving these equations, 
we find the components of the spinor field and metric functions a, Z?, c in terms of volume scale 
T = aba of the BI universe. Finally, in order to find T from Einstein equations and Bianchi identity, 
we deduce three first order ordinary differential equations. Further for convenience we introduce 
a new function v inverse to T, i.e., v = 1/t. 

The fact that the system has the dimension greater than 2, strongly complicates qualitative 
analysis. Note that well known Lorentz system of three ordinary differential equations with poly- 
nomial right hand side with degree less or equal to 2, possesses in some region of parameter space 
chaotic behavior known as a strange attractor and in that region there do not exist first integrals 
(i.e., globally defined invariants). Though the set of singularities is very simple, there exist only 
three singular (fixed) points: two focus and one saddle. The presence of such example does not 
allow us to make an optimistic conclusion on the basis of simple construction of our system (with 
polynomials in the right hand side and absence of singular points the in region of space we are 
interest in, which is even dynamically closed. 

Nevertheless, on the boundary of the the space e = 0, as well as v = (t = +°°), which 
are dynamically closed themselves, the complete classification has been done. The dynamical 
closeness of these planes simultaneously as an obstacle for penetration from positive octant e > 
A V > to the region with negative values. But, there are no singularities, fixed points (there are 
fixed points on the boundary) in the positive octant, we were not able to prove the simplicity of its 
behavior, e.g., presence of first integrals, as well as their absence. 

Thus let us go back to the system (12.11) in details. As it was already mentioned, tt is convenient 
to define a new function v = 1 /t. In this case the obvious singularity that occurs at T = vanishes 
and V = corresponds to T = °o while V = oo to T = 0. The system (12.11) on account of (12.31) takes 
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the form: 



V = -3Hv, (3.1a) 
H = l(35e^//-(l + C)e)-(3//2_£_A)+l(^mv+A(n-2)v""i), (3.1b) 

e = 3//(35e^//-(l + C)e)+4Ae"(3//2-e-A) -4Ae"[mv-Av"j. (3.1c) 



Let us now study the foregoing system of equations in details. 



A. Behavior of the solutions on v = plane 

Let us first study the behavior of the functions H and T on v = plane. The plane v = 
is dynamically invariant, since v|^^q= 0. It should be emphasized that the system in this case 
coincides with one in absence on spinor field and was thoroughly studied in [15]. Nevertheless, 
we write the results obtained in detail. In doing so we rewrite the Eqs. (I3.1bl) and (I3.1cl) in the 
matrix form: 



3H ^ri J \ 3H^-K£-A 

a) By virtu of linear independence of the columns of the matrix of the Eq. (13.21) the critical 
points are the solutions of the equations 

35e^//-(l + C)e = 0, (3.3a) 
3H^-K£-A = 0. (3.3b) 

i.e., they necessarily lie on the parabola (I3.3bl) . In view of 

35 

which follows from (I3.3ah . the Eq. (I3.3bh can be written as 

3fc52ei+2/^ -(l + ^)2e2 + 3A5V^ =0. (3.4) 

The solutions to the system (13.31) will be the roots of the Eq. (13.41) . The quantity of the positive 
roots of Eq. (13.41) according to Cartesian law is equal to the number of changes of sign of the 
coefficients of equations or less than that by an even number. So, for A > and l/2</3<lor 
A < and /3 < 1/2 the number of roots is either 2 or zero. For the remaining cases, i.e., A > 
and /3 > 1 or A > and /3 < 1/2 or A < and /3 > 1/2 there exists only one root. 

In Table 1 classification of qualitatively different types of evolution (phase portrait) depending 
on the parameters /3, A and (1 + Q/B are illustrated. The Fig. a) in Table 1 corresponds to the 
two types, namely /3<l/2or/3 = l/2 and {1 + Q/B < v^3fc, as well as the Fig. i) to the cases 
/3 > 1 or /3 = 1 and {I + Q/B < v^3fc. The Figs, g) and h) of Table 1 cover all the four cases 
for /3 > 1/2. The Figs, d) and e) contain pairs of graphics (case with or 2 singular points; case 
with 1 singular point is also allowable since it possesses the frequency 2, i.e., two singular points 
merges to one saddle-knot and no other qualitative change takes place). They covers 3 cases: d) 
corresponds to /3 < 1 /2; e) corresponds to j8 = 1/2 and {I + Q/B < ^/3k or 1/2 < /3 < 1 and 
j8 = 1/2 and and (1 + Q /B > V3k. 



Anisotropic cosmological models with spinor Geld and viscous fluid ■ 




Table 1. Classification of qualitatively different types of evolution (phase portrait) depending 

the parameters j8, A and {l + Q/B 
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Since, the equation for e only contains 7], the energy density for nontrivial 77 undergoes essen- 
tial changes, whereas H and T remain virtually unchanged. 

The types of critical points lying on the integral curve alternate: . . . saddle, attracting knot, 

saddle So it is sufficient to consider the case with maximum number of roots. Taking into 

account the Eqs. (I2.1cl) and (I3.3bl) let us now calculate 



lim 

£^+°o 3H£ 



lim 



V3B£PVK£+A-e{l + Q 



v^5V/ce(2/5-i) + Ae-2 -{i + Q 



-(1 + C)<0, /3<l/2, 
5v^-(l + C), /3 = l/2, 
+00 >0, /3>l/2. 



(3.5) 



So, the latest critical point for /3 < 1/2 is attracting knot and for /3 > 1/2 is saddle. In case of 
j8 = 1/2 we have either saddle if B\/3k — ( 1 + ^) > and attracting knot otherwise. 

b) It is obvious that if A > the points of intersection of the boundary are the critical points 



H 

£ 



0. 



(3.6a) 
(3.6b) 



c) For H <0 there may exist critical points , if the columns of the matrix of (13.21) are linearly 
dependent. In that case the critical points are the roots of the equation 



and 



3k:(C - l)e + 6k^AB£"+I^ +SkW£^" - 6A = 0, 

2 

H = — kA£". 
3 



In case of T] = the roots of the characteristic equation 

D(H,£) 



D{H, £) 



At 



0, 



are 



(3.7) 
(3.8) 

(3.9) 
(3.10) 



The critical point (H, e) = (0, 2A/[k-(^ - 1)]) is of type divergent focus if A> -9k^^^/[4S{1 + 
Q] or divergent knot if A< -9k^^^ /[4S{1 + Q] . 



1. Integral curves 

For A > the solutions starting from the upper half -plane H > cannot enter into the lower 
one. For A < some of the solutions may enter into the lower half-plane through the segment 
H = and A < < e and never returns back, since H\h=o < 0- 
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B. Behavior of the solutions on £ = plane 

The plane £ = is dynamic invariant, since £ | ^^q= 0. Depending on the sign of H this plane is 
either attractive or repulsive, namely, for // > it is attractive and for // < it is repulsive, since 



-3H{\ + C)<0. 



On £ = plane the system (13.11) takes the form 



V 
H 



-3Hv, 



-3H^ + A+ -{mv + X{n -2)v"-'). 



The system (13.1 II ) has the following integrals: 



3W 
3H^ 



^ 2 . A(n-2)v 
Cv^ + mv+A ^ ^ 



n 



Cv^ + mv + A - ln( v) 



3H^ = Cv^ + mv+A, 



n> 3, 
n = 3, 



(3.11a) 
(3.11b) 



(3.12a) 

(3.12b) 
(3.12c) 



where C is some arbitrary constant. 

The characteristic equation of nontrivial singular points on £ = plane for the system (12.11) 
takes the form 

X{n-2)v"-^ +mv + 2A = 0. (3.13) 

Depending on changes of signs in the sequence of A, m, A it has one, two or no solutions. 

In Table 2 we illustrated the phase-portrait on £ = plane for a positive and a negative A, 
respectively for n = 3. 



A<0 



A>0 




a) 




b) 



Table 2. Classification of qualitatively different types of evolution (phase portrait) on £ = plane 

for n = 3 
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In Table 3 we have graphically illustrated the H — v phase portrait on the e = plane for 
different A. As it was mentioned earlier, here we deal with the multi-parametric system of ordinary 
nonlinear differential equation. In doing so we consider all possible variants independent to their 
physical validity. Therefore, we demonstrate the results obtained for a negative spinor mass (m < 
0). 



A<0 



A = 



A>0 



m < 





a) 



b) 




c) 



m — 




H 0.5 




d) 



e) 



H 0.5 




f) 



m > 



H 0.5 




H 0.5 




g) 



h) 



H 05 




i) 



Table 3. Classification of qualitatively different types of evolution (phase portrait) on e = plane 

for n = 2 



If the right hand side of (13.121) possesses two positive roots with H being positive between 
them, then on the plane e = there occurs closed cycle. It is obvious that there can be no more 
than three roots, hence there cannot be non-concentric cycles. As a result, near the plane e = 
there might be cyclic oscillations. 

The singular point around which the oscillation takes place has H = 0, and therefore, the trajec- 
tory of oscillation partially passes in the region which is attractive to the plane e = and partially 
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in the region that is repulsive. In the long run in the repulsive region at some moment the growth 
of £ becomes dominant. It results in the fact that e becomes infinity within a finite range of time. 



1. Invariants of evolution 

The system (13.11) in absence of viscosity, i.e., under r] —0 and — possesses the following 
first integrals 

Fi = (3.14a) 

F2 = ^ 2 - + ^ (3.14b) 

n-3 

The first of them (|3.14a|) remains to be the first integral even after the introduction of bulk 
viscosity ^. The second one, i.e., Eq. (I3.14bl) under ^ ^ ceases to be the integral of motion. 
Nevertheless, the introduction of bulk viscosity during the course of time generates definite dis- 
placement of the surface given by the formula (|3.14bl) . which allows one qualitatively, i.e., based 
only on the continuity, compile the representation about the possible ways of evolution. 



C. Qualitative analysis of the complete system 

Harnessing the Tables 2, 3 as well as 1, helps one to understand the 3D phase portrait leaning 
on the continuous dependence of the velocity fields of the coordinates v,//, £ of phase space. 

In order to cover the infinite phase space completely, it is mapped on coordinate parallelepiped 
with its axes being the the arc-tangent of the corresponding coordinates. The lower horizontal 
plane always represents the e = plane. 

It should be noted that the introduction of spinor field notably complicates the evolution of the 
system. Contrary to the system in absence of the spinor field, the initial condition with H < 
already does not prevent in many cases thanks to the evolution of volume scale entering the half- 
space H > and thereupon, from the greater value of H repeats the evolution, approaching to the 
V = plane and displaying the classification from the table 1. In the vicinity of the borders £ = 
and V = the integral curves closely repeats the integral curves on the sides, each time at least to 
some extent. 

The general property of all the cases is the fact that in the half-space H > the velocity vectors 
are directed to the £ = plane, while in the other half opposite to it. As a result all he invariant 
curves fall on £ = 0, though not necessarily reach it. 

In the Figs. [T]- [121 we have illustrated the the function inverse proportional to volume scale, 
i.e., v{t) [Fi gs. I1I4I7"TT01| , volume scale T{t) [Figs. I2I5I8I11II and the phase portrait in the v,//,£ 
space [Figs. I3I6I9I12II , corresponding to the positions c,g,h, and i of the Table 2. In doing so we 
used the following values of the parameters: a = 4, j8 = 1, ^ = 0.5, A = 1, B = I and n — 4.The 
positions h and i correspond to the geometrically cyclic regime, but the case h possesses fixed 
point on cyclic integral curve, hence corresponds to the intermediate stage between periodic and 
non-periodic. The positions c and g correspond to the non-periodic evolution. 

The bold black line in the 3D figures [Figs. I3I6I9I12II corresponds to the the functions v(t) 
and T{t) presented in the preceding Fig ures [F igs. 1 11417 II 01215 (STTTll . Figs. [TTIISl and [101111121 
correspond to the case k of Table 1, Fig. l4l5l6l to i and Fig. l7l8l9l to h, respectively. 

The simplest case is illustrated in Figs. IIOII 11121 The integral curves beginning from anywhere 
falls into v = 0. The other peculiarities of this curve are of no significance for the volume scale. 



10 



B. Saha and V. Rikhvitsky 



In case of Fig. I1I2I31 there exists hyperbolic singular point on the £ = plane, so additional to 
the Fig. IIOII lll2l there exists curves tending to the side v = 0°. 

The cyclic trajectories of Figs. I4l5l6l beginning at the vicinity of £ = plane, gradually moving 
away, approach to the £ = 00 plane, and the plane itself in the long run turns out to be attracting 
one with the edges H = ±0°. 

The geometrically cyclic lines in Figs. I7l8l9l in the £ = plane contains fixed point, hence the 
behavior of the corresponding trajectory is very similar to the one given in Figs. IIOII 11121 




FIG. 1: Evolution of func- FIG. 2: Evolution of FIG. 3: 3D view in 

tion inverse to volume volume scale v,//,e space 

scale 




rV. CONCLUSION 

Recently a self consistent system of nonlinear spinor and gravitational fields in the framework 
of Bianchi type-I cosmological model filled with viscous fluid was considered by one of the au- 
thors iliB]- The spinor filed nonlinearity is taken to be some power law of the invariants of bilinear 
spinor forms, namely 7 = 5^ = and J = = (/ij/y^i//)^. Solutions to the corresponding 

equations are given in terms of the volume scale of the BI space-time, i.e., in terms of T = abc, 
with a, b, c being the metric functions. This study generates a multi-parametric system of ordinary 
differential equations Given the richness of the system of equations in this paper a quali- 

tative analysis of the system in question has been thoroughly carried out. A complete qualitative 
classification of the mode of evolution of the universe given by the corresponding dynamic system 
has been illustrated. In doing so we have considered all possible values of the problem parameters 
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FIG. 7: Evolution of func- 
tion inverse to volume 
scale 
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FIG. 10: Evolution of 
function inverse to volume 
scale 




FIG. 8: Evolution of 
volume scale 




D 1 2 3 4 5 



FIG. 11: Evolution 
of volume scale 




FIG. 12: 3D view in 
v,H,e space 



independent to their physical validity and graphically presented the most distinguishable in our 
view results. 
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